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Abstract. Supersymmetry remains compelling theory over 30 years in spite of lack of its discovery. 
It could be already near the corner our days, therefore present and upcoming experiments are crucial 
for constraining or even discovery of the supersymmetry. 

INTRODUCTION 

Since 1970's the Standard Model (SM) based on local gauge invariance principle ex- 
plains all experimental data with a good precision. It is based on SU (3) c x SU {2)i x 
U(l)y non-Abelian Yang-Mills type gauge theory spontaneously broken to SU(3) C x 
U(1)y group. 

Interestingly, at about the same time, when SM has been established, the first ideas 
of Supersymmetry appeared independently around the world: "Extension of the alge- 
bra of Poincare group generators and violation of P invariance" of Golfand and Likht- 
man (1971) id, "Dual theory for free fermions" of Ramon (1971) |0], "Quark model of 
dual pions" of Neveu and Schwarz (1971) |0], "Is the neutrino a Goldstone particle?" 
of Volkov and Akulov (1973) [4]. One should mention separately the work of Wess and 
Zumino, "Supergauge transformations in four-dimensions" (1974) J5|], where the first 
4D supersymmetric quantum field theory has been formulated, which has caused the 
great escalation of the interest of physics community in supersymmetric theories. 

Supersymmetry is the symmetry, which relates bosons and fermions and transforms 
them one into another by acting with supersymmetric fermionic generators Q: 

<2|BOSON) = |FEPvMION), <2|FERMION) = |BOSON) (1) 

One can consider Supersymmetry as one of the most promising attempts to under- 
stand and explain the origin of the fundamental difference between the two classes of 
particles - bosons and fermions. This is already itself quite a reason for the theoretical 
attractiveness of the Supersymmetry. Supersymmetry predicts "mirror" particles to those 
of the SM, which we know. These supersymmetric partners should differ from their SM 
partners by spin 1/2 and have the same mass if the Supersymmetry is unbroken. Since 
it is not the case, and we do not observe this "mirror" supersymmetric world, one can 
conclude that Supersymmetry either does not take place or it is broken. This review is 
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FIGURE 1. The number of pa- 
pers published on Supersymmetry 
versus time 



about the second option, since - as author will try to convince you - there is a big chance 
that Supersymmetry did realized in nature. 

Supersymmetry was invented more then 30 years ago 
but still had not been found yet. However, people are 
publishing thousands of papers on Supersymmetry per 
year till present time and the annual number of the pa- 
pers does not go down as one can see from Fig. 1. What 
makes SUSY being so attractive in spite of the 30 years 
unsuccessful hunting for it?! There are several funda- 
mental reasons for this. 

One of them is that super-Poincare group - the group 
of the Supersymmetry - contains the most fundamental 
set of space-time symmetries as was shown by Haag, 
Lopusanski and Sohnius in 1975 [6]. It contains space- 
time symmetries of the Poincare group and includes 
in addition the supersymmetry transformation, linking 
therefore bosons and fermions. It would be strange if 
the nature did not use this complete set of symmetries. 

Besides this, the requirement for supersymmetric transformation to be local yields 
spin-2 massless gauge filed, the graviton, mediating gravitational interactions 0,11 0- 
Therefore, Supersymmetry could provide unification of all forces in Nature, including 
gravity. 

Supersymmetry also allows to include fermions into Supestring theories which solve 
the problem of non-renormalizability of the gravitational theory and pretend to solve the 
ultimate goal of construction of theory of everything. 

Eventually, the ideas of Supersymmetry - yet a pure theoretical invention - survived 
for more than 30 years because our common belief in unification. 

Miraculous consequence of the Supersymmetry is the prediction of the unification 
of the gauge couplings and solution of the gauge hierarchy problem which are central 
theoretical problems of the SM. In SM one can calculate the evolution of the gauge 
coupling with the energy scale by running respective /3 -functions. Evolution of SU (3), 
SU(2) and U(l) coupling from electroweak (EW) scale up to the 10 15 — 10 16 GeV 
(GUT) scale does not give point-like unification of electromagnetic, electroweak and 
strong forces in case of the SM (Fig. El left frame). In case of minimal Supersymmetric 
extension of the SM (MSSM), SUSY particles modify beta-functions at about 1 TeV 
scale and all three couplings do meet together around Mqut ~ 10 16 GeV (Fig. El right). 



It is worth to stress that attractive features of the supersymmetry mentioned above, 
which solve principal theoretical problems of the SM are only consequences (which has 
been derived many years after the SUSY has been formulated!) of the supersymmetry. 
This fact makes SUSY really attractive from the theoretical point of view. 




SUPERSYMMETRIC GENERALIZATION OF THE STANDARD 

MODEL 

This overview is aimed at giving the current status of SUSY in the light of present 
phenomenological and experimental constraints. Therefore the detailed review of basics 
of Supersymmetry (see e.g. is outside of the scope of this paper 



and will be very brief here. 

Supersymmetry relates bosons and fermions, which transform one into another under 
the action of fermionic generators of supersymmetry algebra Q and Q. Therefore the 
Poincare algebra has to include commutators of 4-momentum and angular momentum 
operators and M^ v with the new fermionic generators Q a and <2« as wei l as me 
anticommutator of Q a and <2cr 

[P^Qa]=0, [M^ v ,Q a ]=-i(an P aQp, {Qa,Qp}=0, {Qa,Qp}=0, {Qa, Qp } = 2(7^ 

(2) 

One can see from © that consequent application of two supersymmetry transformations 
leads to the usual space-time translation generated by P* 1 , which tells us that Supersym- 
metry is a space-time symmetry. Furthermore, making SUSY local, one obtains General 
Relativity, or theory of gravity, or supergravity lZLI&eI]. 

In the simplest case, there is one pair of Q and Q (N = 1 supersymmetry in contrast 
with extended supersymmetries with N > 2). N = 1 is the case of our focus here in 
constructing a supersymmetric version of the Standard Model. 

Superfield formalism is the most convenient way for constructing supersymmetric la- 
grangians. Ordinary fields in this formalism are being extended to superfields which 
are functions of space-time coordinates and anti-commuting Grassmann coordinates 9 a 
and ()P (a,/3 = 1,2). There are two kinds of superfields participating in construction of 
SUSY lagrangian: chiral superfield, 4> and vector superfield, V. Chiral superfield con- 
tains chiral supermultiplet: <I> ~ (A, ty,F), where A(x) is the complex scalar field, 



is the complex Weyl spinor and F(x) is the auxiliary scalar field, which is being elimi- 
nated by equations of motions and which is needed to close the supersymmetric algebra. 
Vector superfield contains V ~ (A,v^,D), where A is a Weyl fermion, is the vector 
(gauge) field and D is the auxiliary field. When constructing SUSY Lagrangian it is con- 
venient to define field strength tensor (an analog of in the non-supersymmetric field 
theory): 

W a = -lD 2 e v D a e- y , and = -JZ) 2 e %r v 

where D and D are covariant derivatives: D a = —d a + i(9o)^ ia d ll , D a = d a — 

and a M = (1,(7'), <7 M = (1,-C7 ! ), are Pauli matrices (i = 1- 3). 

The MSSM Lagrangian should consists of two parts — SUSY generalization of 
the Standard Model and SUSY-breaking part (since we do not observe exact SUSY): 

■^MSSM = -^SUSY + -^breaking, with JZ'sUSY = ^Gauge + ^Gauge —matter "H ^matter — 

£ j( [d 2 TrW a W a + h.c^J + £ f d 2 Qd 2 Q <$>\e Gau z e ^ <D ( + / d 2 8W (3) 

Gauge V"' / Matter-' -* 

where W is a superpotential, which should be invariant under the group of symmetry of 
a particular model and has the following general form: 

J d 2 e + X -mij^j + \y ijk ^^ k ] + h.c. (4) 

The general form of written above superpotential does not forbid violation of baryon 
(B) and lepton (L) numbers. The simultaneous presence of both B- and L- violating 
interactions should be forbidden since it leads to the fast proton decay. To avoid this 
problem one can impose a new discrete symmetry, called i?-parity: R = (— l) 3 ( s_L )+ 25 5 
where S is the spin of the particle. All ordinary Standard Model particles have R = 1, 
the superpartners have R = — 1. R -parity implies that the superparticles can be produced 
only in pairs and that the lightest supersymmetric particle (LSP) is stable, which is a very 
welcome crucial point making LSP the best cold dark matter candidate. In the MSSM 
one assumes that R -parity is conserved. 

When constructing MSSM - minimal supersymmetric extension of the Standard 
Model - one should note that numbers of degrees of freedom of bosons and fermions 
should be the same, and bosonic and fermionic superpartners should have the same 
quantum numbers. Looking at the quantum numbers of the SM particles one concludes 
that new fermions - gauginos (gluino, wino, zino, photino and two higgsinos) - should 
be added to be the superpartners of known bosons and new bosons (scalars) - squarks 
and sleptons - should be added to be the superpartners of known fermions. 

The important feature of the supersymmetric models is the enlarged Higgs sector. 
In particular, Higgs sector of the MSSM has the additional Higgs doublet. One of the 
reasons for it is that in supersymmetric model higgsino has the contribution to the gauge 
anomaly, therefore one needs two higgsinos with opposite hypercharges to realize the 
cancellation of the gauge anomaly. Hence one needs two Higgs doublets with opposite 
hypercharges to make such a cancellation possible. Another reason for the introduction 
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TABLE 1. MSSM particle contents and respective SI/(3) x St/ (2) x 1/(1) quantum numbers 



of additional Higgs doublet is to give masses for both up- and down- type quarks. 
Conjugated Higgs field cannot be used in MSSM (like it was done in SM) since cj>cj>cj>t 
form (which would appear in this analogy) is not chiral superfield and therefore is not 
supersymmetric. 

Table [T] summarizes the particle contents and the quantum numbers of the MSSM. 

The SUSY breaking sector is one of the most essential in MSSM model. We do 
not observe degenerate states of particles and anti-particles, which unavoidably leads 
to conclusion that SUSY should be broken. The implementation of SUSY breaking 
should not spoil the cancellation of the quadratic divergences i.e. SUSY should be 
broken softly. The most general form of the soft SUSY part should include operators 
with mass dimension < 4 iflfjll : 

-&soft= £ m ij A i A J+ E M a (X a X a + l a la) + Y, A ijk^ijkAiAjA k + ^BHiH 2 

scalars gauginos i,j,k 

(5) 

This general form Jtf so f t , which includes many parameters, should be properly con- 
strained in order to address the experimental non-observation of flavor and CP-violating 
processes. 

It is known that none of the fields of the MSSM can develop nonzero vacuum ex- 
pectation value (v.e.v.) to break SUSY without spoiling gauge invariance. In the most 
common scenarios SUSY breaking occur in the hidden sector and propagates to the 
visible sector via messengers. There are several known mechanisms to mediate SUSY 
breaking from hidden to visible sector: gravity mediation (SUGRA), gauge mediation 
(GMSB), anomaly mediation (AMSB) and gaugino mediation (inoMSB). 

In SUGRA scenario lll7Lll8Lll9ll the hidden sector communicates with visible one via 
gravity, leading to the SUSY breaking scale M i ^m^n, where m 3 / 2 is the gravitino 

SUSY 

mass. Scalar masses m,y, gaugino masses M a and trilinear couplings are proportional 
to but can be non-universal in general. In this case one should properly address 
flavor and CP problem. In minimal SUGRA scenario the universality hypothesis of equal 
boundary conditions at the GUT scale greatly reduces the number of parameters down 
to four: m,Q - the common scalar mass, m 1/ / 2 - the common gaugino mass, Aq and B - 



the common tri- and bilinear soft supersymmetry breaking parameters at the GUT scale 
and allows to solve the flavor and CP problem. 

In the case of gauge mediation, messenger is not gravity bu t sup erfield 4> that couples 

l22ll . In this scenario 

(M is the messenger 



to hidden sector and SM fields via gauge interactions [20|, 21, 122I1 

2 ip 

gaugino masses are generated at 1-loop level: M a = -yf^M 
mass scale, while (Fs) is the v.e.v. of singlet scalar superfield of the hidden sector), 

M 



while sfermion masses are generated at 2-loop level: m A 



Remarkable feature of this scenario is that Gravitino is LSP since m.Q ~ 



(Fs) M 



IS 



M Mpi 

suppressed by Plank Mass. Another specific feature of this scenario is that soft masses 
are correlated to the respective gauge couplings. 

msucRA gmsb amsb Anomaly mediation as well as 

Gaugino mediation scenarios are 
based on the paradigm of extra- 
dimensional brane world. Anomaly 
mediated SUSY breaking 11241 12511 is 
generated due to conformal (super- 
Weyl) anomaly at 1-loop level while 
SUGRA mediation effects are ex- 
ponentially suppressed since hidden 
and visible sectors reside on dif- 
ferent branes. In gaugino mediated 




scenario [26, 27], again, hidden and 



FIGURE 3. Sample spectra for SUSY particles for 
SUGRA, GMSB and AMSB scenarios ||H 

visible sectors are on different branes 

while gravity and gauge fields are propagate in the bulk and directly couple fields on 

both branes. As a result, gauginos acquire mass but scalar masses are suppressed and 

may be neglected at the GUT scale. 

The mechanism of SUSY breaking determines the sparticle spectrum and therefore is 

crucial for the SUSY phenomenology and particle searches. Fig. 3 from [23] presents 

the sample spectra for three of four models mentioned above. 



HUNTING FOR SUPERSYMMETRY 

General strategy 

The search for the weak scale Supersymmetry is one of the prime objectives of present 
and future experiments. 

As we discussed above, supersymmetric models can be classified by the mechanism 
for communicating SUSY breaking from the hidden sector to the observable sector. 
Among those scenarios SUGRA could be considered as the most conservative, since it 
requires neither extra dimensions nor new messenger fields while gravity does exist. 

The so-called minimal supergravity (mSUGRA) model has traditionally been the most 
popular choice for phenomenological SUSY analyses. In mSUGRA, it is assumed that 
the MSSM is valid from the weak scale all the way up to the GUT scale Mquj ~ 10 
GeV, where the gauge couplings unify. For this model a simple choice of Kahler metric 



and gauge kinetic function led to universal soft SUSY breaking scalar masses (mo), 
gaugino masses (m l / 2 ) andA-terms (Aq) at Mq UT . This assumption of universality leads 
to the phenomenologically motivated suppression of flavor changing neutral current 
(FCNC) processes. However, there is no known physical principle which gives rise 
to the desired form of Kahler metric and gauge kinetic function, and in general non- 
universal soft breaking mass parameters at the GUT scale are expected II28L 12911 . In 
addition, quantum corrections would lead to large deviations from universality at the 
tree-level Boll . 

We will also require that electroweak symmetry is broken radiatively (REWSB), 
allowing us to fix the magnitude (not the sign) of the superpotential Higgs mass term 
/l to obtain the correct value of Mz- The bilinear soft supersymmetry breaking (SSB) 
parameter B is usually traded for tan/3 (the ratio of Higgs field vacuum expectation 
values). Thus, the parameter set 

mo, mj/2, Ao, tan/3, and sign(fi) (6) 

completely determines the spectrum of supersymmetric matter and Higgs fields. The 
results presented in this review are based on ISASUGRA program, the part of the 
ISAJET ||3l|] package, which calculates the SUSY particle mass spectrum. ISASUGRA 
includes full one-loop radiative corrections to all sparticle masses and Yukawa cou- 
plings, and minimizes the scalar potential using the renormalization group improved 
1-loop effective potential including all tadpole contributions, evaluated at an optimized 
scale choice which accounts for leading two loop terms. Working in the DR regulariza- 
tion scheme, the weak scale values of gauge and third generation Yukawa couplings are 
evolved via 2-loop RGEs to Mquj. At Mquj, universal SSB boundary conditions are 
imposed, and all SSB masses along with gauge and Yukawa couplings are evolved to the 
weak scale M wea ^. Using an optimized scale choice Qsusy = ^ m t L m i R ^ the RG-improved 
one-loop effective potential is minimized and the entire spectrum of SUSY and Higgs 
particles is calculated. There is generally a good agreement between Isajet spectrum and 
publicly available SoftSUSY [32], Spheno 13311 and Suspect [34] codes, as detailed in 
Ref. Q. 

Once the SUSY and Higgs masses and mixings are known, then one should calculate 
observables and compare them against experimental constraints. The most important of 
these are: cold dark matter (CDM) bounds as well as direct and indirect CDM search 
experiments; bounds from various collider experiments (including the future ones) - 
LEP, TEVATRON, LHC, NLC; rare decay processes such as b -> sy, B s -> ji + ji ~, 
fl — > ej\ electric and dipole moments e.g. 8a^ - anomalous magnetic muon moment; 
measurement of the proton life time related to SUSY GUTs theories. 

Cold dark matter bounds. 
Among the most important experiments for SUSY search are those, which confirmed the 
evidence of cold dark matter (CDM) in the Universe and aimed for direct/indirect search 
of CDM candidate. The most direct evidence for CDM in the Universe comes from ob- 
servations of galactic rotation curves. Also binding of galaxies in clusters, matching 
observations of large scale structure with simulations, gravitational microlensing, bary- 
onic density of the Universe as determined by Big Bang nucleosynthesis, observations 
of supernovae in distant galaxies, and measurements of anisotropics in the cosmic mi- 



crowave background radiation (CMB) can be considered as a confident confirmation of 
CDM (for reviews see e.g. jl36U37lO . In particular, the Wilkinson Microwave Anisotropy 
Probe (WMAP) [13 8L l39fl collaboration has extracted a variety of cosmological parame- 
ters from fits to precision measurements of the CMB radiation. It was found that 



a 



total 



1.02 ±0.02 



(7) 



where Q. total = p to tal/pc, Pwtai is the matter-energy density of the Universe, p c (G N is 
the Newton constant) the critical density of the closure of the Universe, defined within 
Friedman-Robert- Walker (FRW) framework. H is the Hubble parameter usually defined 
through the rescaled Hubble parameter h as H — 100/z Km/(s Mpc). So, with a good 
precision, Universe has been measured to be flat. The properties of the Universe are 
characterized by the density of baryons matter density (Q. m ), vacuum energy (^a) 
and the expansion rate (h) which are measured to be: 



^ = 0.044 ±0.004, a m = 0.27 ±0.04, Q A = 0.73 ±0.04, /i = 0.71 
From the WMAP results, one derives the following value for CDM: 



+0.04 
-0.03- 



(8) 



Q-CDMh 



^ h h 2 = 0A\26t S(t° S) at68(95)%CL. (9) 



The energy content of the Universe according to WMAP data is schematically presented 
in Fig. 4. There exists a number of hypothetical candidate elementary particles to fill the 
role of CDM. A particularly attractive candidate for CDM is the lightest neutralino in 
7?-parity conserving supersymmetric models EoL l4lh . 

■ ^^^^^^^M The evolution of the number density of 

^B^^^^ _^^^y^^5* j^^BB supersymmetric relics in the universe is de- 
WF BjCtf *B scribed by the Boltzmann equation as formu- 

¥ H lated for a FRW Universe. For calculations 

I ^^T^l including many particle species, such as the 

L ■ ■TiQy^vBlMJ case wnere co-annihilations are important. 

Wk liTifivj^,^ '-i J there is a Boltzmann equation for each particle 

^^k\ species H42I1 . the equations can be combined to 
| obtain 



FIGURE 4. The energy content of the Uni- 
verse according to WMAP data (source: 
http://map.gsfc.nasa.gov/! 



— = -3Hn - (o ef fv) [n - n eq ) 



(10) 



2 j 

where n = YJtL\ n U n eq,i = K2 (y-) and the sum extends over the /V particle species 
contributing to the relic density, with n ; being the number density of the i-th species, Kj 
is a modified Bessel function of the second kind of order j. The quantity (o e ffv) is the 
thermally averaged cross section times velocity which is the key point in calculating 
relic density including co-annihilation processes with Cih 2 \ l\ ^ F ( o P ffv )dx\, where 
one integrates (a e ffv) from to freeze-out temperature, Xf II43L l42l l44fl . The main 
challenge in CDM calculation is the evaluation of all possible channels for neutralino 
annihilation to SM and/or Higgs particles, as well as all co-annihilation reactions with 
relativistic thermal averaging. For our calculation of the neutralino relic density, we use 
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FIGURE 5. Regions of neutralino relic density in the mo vsM]^ plane for Aq = and tan/3 = 30(55). 
Generic regions which are in agreement with WMAP results: 1) bulk annihilation region, 2) stau co- 
annihilation region, 3) hyperbolic branch/focus point (HB/FP) region and 4) A-annihilation funnel region 
(only right plot for high tan/3). 



the IsaReD program interfaced with Isajet. There are also public programs for CDM 
calculations such DarkSUSY S H3] and micrOMEGAs H HH (as well as many 
private ones) which are in a good agreement with IsaReD. In most of the parameter space 
of the mSUGRA model, it turns out that a value of ^ 2 we U beyond the WMAP bound 
is generated. Only certain regions of the mSUGRA model parameter space give rise to a 
relatively low value of ^z^ 2 m accor d with WMAP measurements. Fig.|5]presents such 
generic regions in mo versus m\/2 plane for \x > 0, Aq = tan/3 = 30,55 (left, right). 
The dark (red) shaded regions are excluded by theoretical constraints (lack of REWSB 
on the right, a charged LSP in the upper left). The blue (dark-gray) shaded region 
(0. 129 > Q.j h 2 > 1) and, furthermore, unshaded region (% h 2 > 1) are excluded, since 
their too high relic density would lead to overclosure of the Universe contradicting with 
WMAP data. The green (gray) shaded regions are in exact agreement with WMAP 
constraints (Eq.|3) and could be classified in several classes: 

1. The bulk annihilation region at low values of mo and m 1/ / 2 , where neutralino pair 
annihilation occurs at a large rate via ^-channel slepton exchange. 

2. The stau co- annihilation region at low mo where ~ so that Zis may co- 
annihilate with fis in the early universe ll50ll5lll . 

3. The hyperbolic branch/focus point (HB/FP) region [52i 23 > 24 > 22] at large mo near 



the boundary of the REWSB excluded region where \fi\ becomes small, and the 
neutralinos have a significant higgsino component, which facilitates annihilations 
to WW and ZZ pairs. 

4. The A-annihilation funnel, which occurs at very large tan/3 ~ 45 — 60 ll56l 22, 



l58l l59t l60l l61h . In this case, the value of m^ ~ 2m= . An exact equality of the 



mass relation isn't necessary, since the A width can be quite large (Fa ~ 10 — 50 
GeV); then 2/n^ can be several widths away from resonance, and still achieve a 

large Z\Z\ — > A — > // annihilation cross section. The heavy scalar Higgs H also 
contributes to the annihilation cross section. 

In addition, there exists a region of neutralino top-squark co-annihilation (for very 
particular Aq values) and a light Higgs h annihilation funnel (at low m\ii values). In the 

light-gray (yellow) region neutralino relic density is too low (Q^h 2 < 0.096) to give a 
good fit for WMAP data. But, strictly speaking, this region is not excluded since there 
could be additional particles (e.g. axions or particles from various models with extra- 
dimensions) contributing to CDM relic density. Therefore, hereafter we treat green and 
yellow regions of Fig.|5]as one, allowed ("green") region. 

Constraints from LEP2 searches 
The LEP2 collaborations have finished taking data and collected the rich statistics for 
e + e~ CM energies ranging up to y/s ~ 208 GeV. Searches for superpartners at LEP2 
gave negative results and led to the following constraints: 



m% > 103.5 GeV B62J], m- e > 99 GeV provided m ? -m^ < 10 GeV (63J (11) 



The LEP2 experiments also searched for the SM Higgs boson. In addition to finding 
several compelling signal candidates consistent with ~ 115 GeV, they set a limit 
niH SM > 114.1 GeV 16411 . In our mSUGRA parameter space scans, the lightest SUSY 
Higgs boson h is almost always SM-like. The exception occurs when the value of 
becomes low, less than 100— 150 GeV. This Higgs mass LEP2 limit requires radiative 
corrections to the light Higgs mass (which is below Z-boson mass at the tree level) 
to be quite large. Those radiative corrections - 8Mh - are logarithmically sensitive to 
the SUSY scale (M SUSY ) and this LEP2 limit pushes this to level of ~ 1 TeV. At the 
same time 8Mj 2 mf op . In the light of recent m top measurements from DO ll65ll which 
shifted up top-quark mass to be 178.0±4.3 GeV (as a world averaged value) led to some 
relaxation of the lower limit on Msusy- 

The b — > sy branching fraction 
The branching fraction BF{b — > sy) has recently been measured by the BELLE ll66ll . 
CLEO0 

and ALEPH [68] collaborations. Combining statistical and systematic errors 
in quadrature, these measurements give (3.36 ±0.67) x 10~ 4 (BELLE), (3.21 ±0.43) x 
10~ 4 (CLEO) and (3. 1 1 ± 1 .07) x 10~ 4 (ALEPH). A weighted averaging of these results 
yields BF(b -> sy) = (3.25 ±0.37) x 10~ 4 . The 95% CL range corresponds to ±2<7 
away from the mean. To this we should add uncertainty in the theoretical evaluation, 
which within the SM dominantly comes from the scale uncertainty, and is about 10%. 1 
Together, these imply the bounds, 

2.2 x 10~ 4 <BF(b^sy) < 4.33 x 10~ 4 (12) 



We caution the reader that the SUSY contribution may have a larger theoretical uncertainty, particularly 
if tanj3 is large. 



Other computations of the range of BF(b — > sy) include for instance Ellis et al. [|6S 
2.33 x 10- 4 <BF(b^sy) <4.15x 10 4 ,andDjouadiefa/. & 2.0 x 10~ 4 <BF(b^ 
sy) < 5.0 x 10~ 4 . In our study, we simply show contours of BF(b — > sy) of 2, 3, 4 and 
5 x 10 -4 , allowing the reader the flexibility of their own interpretation. 

The calculation of BF(b — > sy) used here is based upon the program ofRef. ED, 



7211 . That calculation uses an effective field theory approach to evaluating radiative 
corrections to the b — > sy decay rate. In our calculations, we implement the running 
^-quark mass including SUSY threshold corrections as calculated in ISASUGRA; these 
effects can be important at large values of the parameter tan /3 17 3l 17411 . Once the relevant 
operators and Wilson coefficients are known at Q = My/, then the SM WCs are evolved 
down to Q = mt via NLO RG running. At m^, the BF(b — > sy) is evaluated at NLO, 
including bremsstrahlung effects. Our value of the SM b — > sy branching fraction yields 
3.4 x 10~ 4 , with a scale uncertainty of 10%. 
Muon anomalous magnetic moment 

The muon anomalous magnetic moment a u = ^ 2 2 " >M nas been recently measured to 
high precision by the E821 experiment Il75ll for the negative muon along with earlier 
results on the positive muon [76]. In addition, theoretical determinations of (g — 2)^ 
have been presented by Davier et al. ¥Fh and Hagiwara et al. iFsIl which use recent 
data on — > hadrons at low energy to determine the hadronic vacuum polarization 
contribution to the muon magnetic moment. Combining the latest experiment and theory 
numbers, we find the deviation of a^ to be: 

Aa^ = (27.1 ±9.4) x 10" 10 (Davier et al.) (13) 
Aa^ = (24.7 ± 9.0) x 10" 10 (Hagiwara et al.). (14) 

The Davier et al. group also presents a number using x decay data to determine 
the hadronic vacuum polarization, which gives Aa^ = (12.4 ±8.3) x 10 -10 , i.e. nearly 
consistent with the SM prediction. However, there seems to be growing consensus 
that the numbers using the e + e data are to be trusted more, since they offer a direct 
determination of the hadronic vacuum polarization. The ~ 3a deviation in a^ using the 
e + e~ data can be explained in a supersymmetric context if second generation sleptons 
(smuons and muon sneutrinos) and charginos and neutralinos are relatively light. 

#s — > decay 

While all SUSY models contain two doublets of Higgs superfields, there are no tree- 
level flavor changing neutral currents because one doublet H u couples only to T3 = 1/2 
fermions, while the other doublet Hj couples just to T3 = — 1/2 fermions. At one loop, 
however, couplings of H u to down type fermions are induced. These induced couplings 
grow with tan /3 . As a result, down quark Yukawa interactions and down type quark mass 
matrices are no longer diagonalized by the same transformation, and flavor violating 
couplings of neutral Higgs scalars h, H and A emerge. Of course, in the limit of large 
niA, the Higgs sector becomes equivalent to the SM Higgs sector with the light Higgs 
boson h = Hsm, and the flavor violation decouples. The interesting thing is that while 
this decoupling occurs as tha — > 00 , there is no decoupling for sparticle masses becoming 
large. 



An important consequence of this coupling is the possibility of the decay B s — > , 
whose branching fraction has been experimentally bounded by CDF [79] to be: 

BF(B S < 2.6 x 1(T 6 , (15) 

mediated by the neutral states in the Higgs sector of supersymmetric models. While 
this branching fraction is very small within the SM (BFsm(B s — * M + M ) — 3.4 x 10~ 9 ), 
the amplitude for the Higgs-mediated decay of B s grows as tan 3 /3 within the SUSY 
framework, and hence can completely dominate the SM contribution if tan /3 is large. 
Several groups [80, 8jJ, lS2J] and recently 118311 have analyzed the implications of this 
decay within the mSUGRA framework. We also present improved b — > sy branching 
fraction predictions in accord with Il83ll for the current ISAJET release. This constraint 
is important only for very large values of tan/3. 

mSUGRA constraints. 
Within the mSUGRA framework, the parameters mq and n are the most important 
for fixing the scale of sparticle masses. The mo — m\/2 plane (for fixed values of other 
parameters) is convenient for a simultaneous display of these constraints, and hence, 
of parameter regions in accord with all experimental data. Physicists interested in the 
mSUGRA model may wish to focus their attention on these regions. 

Our results for combined constraints mentioned above are presented in Fig. |5] for 
mSUGRA model in the mO vs m\n plane for Ao = and tan/3 =30 (55) in the left 
(right) frame, respectively. Top raw of frames presents constraints themselves, while the 
bottom raw presents corresponding % 2 pattern for frames above. % 2 is formed from Aa^, 
BF(b -> sy) and Sl^h 2 , as given in Eqn's (1911 211 41) . 

Allowed CDM relic density regions are the green shaded ones, LEP2 excluded regions 
denoted by yellow color. Black contours denotes 8a ^ values (30, 10, 5, 1 x 10~ 10 - from 
bottom to top) and blue contours denote BF(b — > sy) values (2, 3 x 10~ 4 - from bottom 
to top). In the left upper frame of Fig. |5] one can observe the tension between 8a n and 
BF(B — > sy). The point is that while BF(b — ► sy) favors large third generation squark 
masses to suppress SUSY contributions to b — > sy decay, Aa^ experimental value favors 
relatively light second generation slepton masses, to give a significant (g — 2) „ deviation 
from the SM value. This is clearly reflected in the respective (left bottom) frame where 
the % 2 quantity is presented. 

In case of low and intermediate values of tan/3 almost all the mo vs. m\/2 plane 
has very large % 2 . This arises because in general an overabundance of dark matter is 
produced in the early universe, and the relic density ^■j.b 1 is beyond WMAP limits. 
There is a very narrow sliver of yellow at m^n ~ 150 GeV (just beyond the LEP2 
limit) where 2m^ ~ m^, and neutralinos can annihilate through the narrow light Higgs 
resonance. In addition, there is an orange/yellow region at high mq at the edge of 
parameter space (the HB/FP region), with an intermediate value of % 2 . In an earlier 
study [84], this region was found to have a low x 1 value. In the present situation, 
however, the 3o deviation from the SM of a^ tends to disfavor the HB/FP region. In the 
HB/FP region, sleptons are so heavy (typically 3-5 TeV), that SUSY contributions to a^ 
are tiny, and the prediction is that a^ should be in near accord with the SM calculation. 
The remaining green region is the narrow sliver that constitutes the stau co-annihilation 
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FIGURE 6. Constraints for mSUGRA model in the ;«o vs mi/2 plane for Aq — and tan/3 = 30(55) 
(top raw). Top raw: allowed CDM relic density regions are the green shaded ones, LEP2 excluded regions 
denoted by yellow color. Black contours denotes Sa^ values (30, 10,5, 1 x 10~ 10 - from bottom to top) 
and blue contours denote BF(b — > sy) values (2,3 x 10~ 4 - from bottom to top). Magenta contours for 
BF(Bs — > are relevant only for very high tan/3 = 55 (7.5, 1 x 10~ 7 - from bottom to top). The 

bottom raw of figures presents \p/} for the respective parameters of the mSUGRA model. The green 
regions have low x 2 , while red regions have high % 2 . Yellow is intermediate. 



region, barely visible at the left hand edge of parameter space adjacent to where fi 
becomes the LSP. 

Once we move to very large tan /3 values, as shown in the right set of frames of Fig.|6| 
then the A-annihilation funnel becomes visible, and some large regions of moderately 
low X 2 appear around mo, mj/2 ~ 500, 600 GeV and also at 1500, 200 GeV. While 



the A-annihilation funnel extends over a broad region of parameter space, the upper and 
lower ends of the funnel are disfavored: basically, if sparticles become too heavy (the 
upper end), then Aa^ becomes too small, while if sparticles become too light (the lower 
end), then BF(b — > sy) deviates too much from its central value. 

One comes to conclusions II85I1 that for the mSUGRA model almost all of parameter 
space is excluded or at least disfavored by the combination of the WMAP Q.~ h 2 limit, 

the new Aa^ value, and the BF(b — > sy) value. The f^/z 2 constraint only allows 
the several regions of parameter space mentioned above, while BF{b — > sy) favors 
large third generation squark masses to suppress SUSY contributions to b — > sy decay, 
and Aa^ favors relatively light second generation slepton masses, to give a significant 
deviation of (g — 2)u from the SM value. The only surviving regions with relatively 

low \p)ft- ~ 2 are the stau co-annihilation region, and intermediate portions of the A- 
annihilation funnel at very large values of tan /3 . One should note, however, that if the 
hadronic vacuum polarization determination using x decay data turn out to be correct, 
then the HB/FP region will appear in a more favorable light! 

Reach of LHC and LC in Dark Matter Allowed Regions of the mSUGRA Model 
The reach of the CERN LHC for SUSY in the mSUGRA model has been calculated in 
Ref. [86] assuming 100 fb 1 of integrated luminosity. Briefly, sparticle pair production 
events were generated for many mSUGRA model parameter choices in the niQ vs. m^n 
plane for various tan/3 values. A fast LHC detector simulation (CMSJET) was used, and 
cuts were imposed to extract signal rates in a variety of multilepton plus multijet plus 
missing transverse energy channels. Backgrounds were also calculated from a variety of 
QCD and vector boson production processes. A large set of selection cuts were used to 
give some optimization over broad regions of parameter space. It was required to have 
at least a 5o signal over background, with at least 10 signal events in the sample. 

The reach of the CERN LHC is shown in Fig. for the case of tan j8 = 30 (55) left 
(right), /i > 0, Ao = and m t = 175 GeV. The dark shaded (red) regions are disallowed 
by lack of radiative electroweak symmetry breaking (REWSB) (right hand side) or 
presence of a stau LSP (left hand side). The light gray (yellow) region is excluded 
by LEP2 chargino searches (mti > 103.5 GeV), while the region below the yellow 
contour gives m/ 7 < 1 14.4 GeV, in contradiction of LEP2 SM Higgs searches (here, 
the SUSY h Higgs boson is essentially SM-like). The medium gray (green) regions 
have CLcDMh 2 < 0.129, and are allowed by WMAP. The broad HB/FP region is seen 
on the right-hand side, while the stau co-annihilation region is shown on the left- 
hand side. At the edge of the LEP2 excluded region is the light Higgs annihilation 
corridor. The reach of the Fermilab Tevatron via the trilepton channel is also shown ll87Tl . 
assuming a 5c signal over background for 10 fb 1 . The reach of the CERN LHC for 100 
fb 1 of integrated luminosity is shown by the contour labeled "LHC". It extends from 
m i/2 ~ 1400 GeV (corresponding to a value of nig ~ 3 TeV) on the left-hand side, to 
m i/2 ~ ^00 GeV (corresponding to mg ~ 1.8 TeV) on the right-hand side. In particular, 
for those values of tan/3, the LHC reach covers the entire stau co-annihilation region, 
plus the low nil j2 portion of the HB/FP region. The outer limit of the reach contour is 

mainly determined by events in the Ef lss + jets channel, which arises from gluino and 
squark pair production, followed by hadronic cascade decays. 
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FIGURE 7. Parameter space of mSUGRA model for tan/3 = 30(55)left(right), A = and ju > 0, 
showing the reach of the Fermilab Tevatron, the CERN LHC and a 0.5 and 1 TeV linear e + e~ collider for 
supersymmetry discovery. 



We also show in the plot the reach of a y/s = 500 and 1000 GeV LC, assuming 
100 fb 1 of integrated luminosity HSSfl . Events were generated using Isajet 7.69, and 
compared against various SM backgrounds. The left-most portion of the reach contour 
arises where selectron and smuon pair production are visible, while the main portion (flat 
with m 1/ / 2 ) arises due to chargino pair searches. An additional reach is gained between 

these two regions by searching for e + e~ — > X2X1 production, followed by ~^ l\bb 
decay. In addition, in Ref. [ 89], additional reach can be gained by searching for stau pair 
events, although two photon backgrounds must be accounted for, due to the low energy 
release in the stau co-annihilation region. 

While a 500 GeV LC can cover only a portion of the stau co-annihilation region, a 1 
TeV LC can cover the entire region, at least for this value of tan/3. As one moves into 
the HB/FP region, the LC retains a significant reach for SUSY, which in fact extends 
beyond that of the CERN LHC! It is significant that this additional reach occurs in a DM 
allowed region of parameter space. In the HB/FP region, the superpotential \x parameter 
becomes small, and the lightest chargino and neutralino become increasingly light, with 
increased higgsino content. In fact, the decreasing mass gap between %\ an d %? makes 
chargino pair searches difficult at a LC using conventional cuts because there is so little 
visible energy release from the chargino decays. In Ref. [88, 89], we advocated cuts that 
pick out low energy release signal events from SM background, and allow a LC reach 
for chargino pairs essentially up to the kinematic limit for their production. In this case, 
it is important to fully account for 77 — > ff backgrounds, where / is a SM fermion. 

The case study of Ref. [88] also examined the HB/FP region with parameters mo = 
2500 GeV, m 1/2 = 300 GeV, A = 0, tan/3 = 30, jl > and m, = 175 GeV, i.e. in the 
HB/FP region. In this case, chargino pair events were selected from the l£+jets +Ejf ss 
channel, and the diiet mass distribution was used to extract the value of m v + and m^o 
at the 10% level. The mass resolution is somewhat worse than in previous case studies 



in the literature because the charginos undergo three-body rather than two-body decays, 
and no sharp edges in energy distributions are possible. Nonetheless, the measured value 
of chargino and neutralino mass, along with a measure of the total chargino pair cross 
section, was enough to determine the SUSY parameters M2 and /i to 10-20% precision. 
The results, shown in Fig. 8, demonstrate that /1 < Mi, which points to a X\ an d X\ 
which are higgsino/gaugino mixtures, as is characteristic of the HB/FP region. 

In conclusion to this section, one should 
stress that the CERN LHC and an e+e LC 
are highly complementary to each other in 
exploring the dark matter allowed param- 
eter space of the mSUGRA model. LHC 
covers the stau co-annihilation region (com- 
pletely for tan/3 < 40) as well as the H, A 
funnel region (much of which is typically 
beyond the maximum reach of a LC). How- 
ever only the lower part of the HB/FP re- 
gion can be covered by the LHC. On the 
other hand, as we have demonstrated, LCs 
can probe much of the upper part of the 
HB/FP region with the new proposed cuts. 
Therefore, the combination of the LHC and 
a TeV scale LC can cover almost the entire 
parameter space of the mSUGRA scenario. 

Direct and Indirect CDM searches 
Besides collider experiments there exist both direct and indirect non-accelerator dark 
matter search experiments, which are ongoing and proposed. Those experiments as will 
be discussed below play crucial role in restricting SUSY parameter space. 
Direct CDM search experiments are aimed to look at neutral dark matter candidates 
scattering off nuclei. Direct dark matter detection has been recently examined by many 
authors 19011 . and observable signal rates are generally found in either the bulk annihi- 
lation region, or in the HB/FP region, while direct detection of DM seems unlikely in 
the A-funnel or in the stau co-annihilation region. Early limits on the spin-independent 
neutralino-nucleon cross-section (05/) have been obtained by the CDMS 119 ill . EDEL- 
WEISS l&and ZEPLIN1 Q groups, while a signal was claimed by the DAMA col- 
laboration Il94ll . Collectively, we will refer to the reach from these groups as the "Stage 
1" dark matter search. Depending on the neutralino mass, the combined limit on Osi 
varies from 10 5 to 10 pb. This cross section range is beyond the predicted lev- 
els from most supersymmetric models. However, experiments in the near future like 
CDMS2, CRESST2 B, ZEPLIN2 and EDELWEISS 2 (Stage 2 detectors) should have 
a reach of the order of 10~ 8 pb. In fact, the first results from CDMS2 have recently 
appeared, and yield a considerable improvement over the above mentioned Stage 1 re- 
sults S. Finally, a number of experiments such as GENIUS 0, ZEPLIN4 S and 
XENON [99] are in the planning stage. We refer to these as Stage 3 detectors, which 
promise impressive limits of the order of Osi < 10~ 9 - 10~ 10 pb, and would allow the 
exploration of a considerable part of parameter space of many supersymmetric models. 
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FIGURE 8. Determination of SUSY parameters 
from examining chargino pair production at a 
y/s = 0.5 TeV LC, for the HB/FP mSUGRA point 
listed in the text. 



In particular, the Stage 3 direct DM detectors should be able to probe almost the entire 
HB/FP region of mSUGRA model parameter space. We note here in addition that the 
Warm Argon Program (WARP) [ 100] promotes a goal of detecting neutralino-nucleus 
scattering cross sections as low 10~ u pb. 



Indi rect detection of neutralino dark matter may occur via (for review, see e.g. [ 101L I102 , 

floi 

1 . observation of high energy neutrinos originating from Z\Z\ annihilations in the core 
of the sun or earth 

2. observation of y-rays originating from neutralino annihilation in the galactic core 
or halo and 

3. observation of positrons or anti-protons originating from neutralino annihilation in 
the galactic halo. 

The latter signals would typically be non-directional due to the influence of galactic 
magnetic fields, unless the neutralino annihilations occur relatively close to earth in 
regions of clumpy dark matter. The indirect signals for SUSY dark matter have been 
investigated in a large nu mber of papers, and computer codes which yield the various 
signal rates are available lll04ll47ll . Recent works find that the various indirect signals 
occur at large rates in the now disfa vored bulk annihilation region, and also in the 
HB/FP region ill05Lll06l1 . In Ref. [107], it was pointed out that the A annihilation funnel 
can give rise to large rates for cosmic ys, e + s and ps. However, neutralino-nucleon 
scattering cross sections are low in the A annihilation funnel, so that no signal is expected 
at neutrino telescopes, which depend more on the neutralino-nucleus scattering cross 
section than on the neutralino annihilation rates. 

There is a bunch of various experiments for indirect dark matter detection. Neutrino 
telescopes such as Antares or IceCube are aimed to detect muon neutrinos from neu- 
tralino annihilation in the core of the earth or sun via — ► /! conversions. The Antares 
v telescope should be sensitive to E jr > 10 GeV; it is in the process of deployment 
and is expected to turn on in 2006 [i08]. It should attain a sensitivity of 100— 1000 
/is/km 2 /yr. The IceCube v telescope is also in the process of deployment at the south 
pole fl 1 OOL 1 1 1 dl . It should be sensitive to E jl > 25 — 50 GeV, and is expected to attain a 
sensitivity of 10 — 100 /is/km 2 /yr. Full deployment of all detector elements is expected 
to be completed by 2010. 

Neutralino could also annihilate in the galactic halo giving rise gamma rays, 
positrons or anti-protons fluxes. The y rays can be detected down to sub-GeV energies 



with space-based detectors such as EGRET 111 1111 or GLAST [112]. Ground based 
arrays require much higher photon energy thresholds of order 20— 100 GeV. Experi- 
ments such as GLAST should be sensitive to rates of order 10~ 10 ys/cm 2 /sec assuming 
Ey > 1 GeV. In fact, it has recently been suggested that the extra-galactic gamma ray 
background radia tion as m easu red by EGRET is well fit by a model of neutralino 



annihilation 1 113 . 114 1. In llll5ll detailed study has been done explaining not only the 



rate of SUSY contribution to the EGRET data but also the shape of various distribution 
through the non-trivial shape of CDM halo. 

Positrons would arise as decay products of heavy quarks, leptons and gauge bosons 
produced in neutralino annihilations. Space based anti-matter detectors such as 
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FIGURE 9. A plot of the reach of direct, indirect and collider searches for neutralino dark matter in the 
mo vs. m\/2 plane, for Aq = 0, tan/3 = 30(55)left(right) and jj. > 



Pamela llllfjll and AMS-02 1 1 17 1 will be able to search for anomalous positron 
produc tion f rom dark matter annihilation. The cosmic positron excess as measured by 
HEAT 1 1 181 has been suggested as having a s ource in ga lactic halo neutralino annihila- 
tions ll 1 1 0L 1 1 14f] . It is suggested by Feng et al. lll05lll06|l that a reasonable observability 
criteria is that signal-to-background (S/B) rates shou ld be gre ater than the 1-2% level. 



To calculate the S/B rates, we adopt fit C from Ref. lll05Lll06ll for the E 2 d® e + jdQ.dE 
background rate: 

E 2 d® e+ /dndE = 1.6 x 1(T 3 £~ L23 , (16) 

where E is in GeV. 

Anti-protons may also be produced in the debris of neutralino annihilations in th e 
galactic halo. Such anti-protons have been measured by the BESS collaboration ill2(lll . 
The differential flux of anti-protons from the galactic halo, d&p/dEpdQ., as measured 
by BESS, has a peak in the kinetic energy distribution at Ep ~ 1.76 GeV. The height of 
the peak at Ep ~ 1.76 GeV is ~ 2 x 10~ 6 p/GeV/cm 2 /s/sr. Signal rates in the range of 
10~ 7 — 10~ 6 p/GeV/cm 2 /s/sr might thus provide a benchmark for observability. 

Our results for direct/indirect DM search experimental reach for the mSUGRA model 
are shown in Fig. [9] in the mo vs. m\ji plane for tan/3 = 30 and 55 (left and right 
respectively), Ao = and \x > 0. If addition to collider reach presented in the previous 
figures, here w e show contours of 

• Stage 3 direct detection experiments (05/ > 10~ 9 pb; black contour), 

• reach of IceCube V telescope with <J>™"(/i) = 40 jUs/km 3 /yr and E^ > 25 GeV 
(magenta contour), 



• the 4>(y) = 10 10 ys/cm 2 /s contour with Ey > 1 GeV in a cone of 0.001 sr directed 
at the galactic center (dark blue contour), 

• the S/B > 0.01 contour for halo produced positrons (blue-green contour) and 

• the anti-proton flux rate <&(p) = 3 x 10~ 7 ps/cm 2 /s/sr (lavender contour). 

As noted by Feng et al. lll05l Il06ll . all these indirect signals are visible inside some 
portion of the HB/FP region, while none are visible in generic DM disallowed regions 
(under the assumed smooth halo profiles). The intriguing point is that almost the entire 
HB/FP region can be explored by the cubic km scale IceCube v telescope! It can also 
be explored (apparently at later times) by the Stage 3 direct DM detectors. Taking into 
account relative time scales of the various search experiments, if SUSY lies within the 
upper HB/FP region, then it could be discovered first by IceCube (and possibly Antares), 
with a signal being later confirmed by direct DM detection and possibly by the TeV scale 
linear e + e~ collider. There is also some chance to obtain indirect y, e + and p signals in 
this region. Notice that if instead SUSY lies within the stau co-annihilation corridor, 
then it will be discovered by the LHC, but all indirect detection experiments will find 
null results in their DM searches. 



BEYOND MSUGRA: NORMAL MASS HIERARCHY 

We see that Aa^ favors light second-generation sleptons, while BF[b — ► sy) prefers 
heavy third generation squarks. Since this situation is hard to realize in the mSUGRA 
model, this could be an indication that one must move beyond the assumption of univer- 
sality, wherein each generation has a common mass at Q = Mquj. Therefore scenario 
in which first and second generation scalars remain degenerate, while allowing for a 
significant splitting with third generation scalars is well motivated. In this case, heavy 
(multi-TeV) third generation scalars are preferred by BF(b — > sy) constraints, while 
rather light first and second generation scalars are preferred by Aa^. The scenario is 
called the normal scalar mass hierarchy (NMH). The parameter set of the mSUGRA 
model is expanded to the following values: H85I1 

m (l), m (3), m H , m 1/2 , A , tan/3, signal). (17) 

where mo(l) is the common scalar mass of all first and second generation scalars at 
Q = Mquj, while mo(3) is the common mass of all third generation scalars at Mquj. 
The above parameter set is well motivated in 5O(10) SUSY GUT models, where the two 
MSSM Higgs doublets typically occupy a 10 of 5O(10), and each generation of scalars, 
along with a SM gauge singlet /V occupies the 16 dimensional spinorial representation of 
5O(10). The step of breaking generational universality must be taken with some caution, 
since in general it can lead to violations of constraints from FCNC processes. Splitting 
the third generation from the first and second can potentially lead to violations of FCNC 
processes. One of the main experimentally measured bounds on FCNC processes in this 
case comes from B Q H —B Q L mass splitting. One can show that Amg is much less restrictive 
than the kaon case, for both low and high squark masses. Moreover, even if squark mass 
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FIGURE 10. Left frame: ^fx} values in the m (3) vs. m 1/2 plane for m (l) = 100 GeV, tan/3 = 30, 
Ao = and ji > in NMH SUGRA scenario. The green regions have low x 2 , while red regions have 
high x 2 . Yellow is intermediate. Right frame: constraints for NMH SUGRA model in the mo(3) vs. 
m\ji plane for mo(l) = lOOGeV, Aq = and tan/3 = 30. Allowed CDM relic density regions are the 
green shaded ones, LEP2 excluded regions denoted by blue color. Black contours denotes da^ values 
(30, 10,5, 1 x 10~ 10 - from bottom to top) and blue contours denote BF(b — > sf) values (2,3 x 10~ 4 
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splittings are very large at the scale Q = Mquj, the weak scale mass splittings are much 
smaller and practically whole parameter space of NMH SUGRA is allowed [85]. 

Scan in the SUGRA parameter space of Eq. [T7] and respective x 2 analysis show 
the following preferences II85I1 : m# ~ mo(3), mo(l) C mo(3), mo(l) ~ — 400 GeV, 
mo (3) ~ 500 — 3000 and large iuq ( 1 ) — mo (3) mass splitting and that mainly large tan /3 
is preferred if the generational mass splitting is small (which takes us back towards 
the mSUGRA case). If may also be pointed out that SUSY IMH models are greatly 
disfavored. The above scan motivates to reduce the number of parameters and choose 
rriH = mo (3) which brings us to the minimal extension of mSUGRA space by just 
one additional parameter mo — > [mO(l),mo(3)]. In Fig. flOT left) we present the \py} 
values in the mo(3) vs. ni\ii plane for mo(l) = 100 GeV, tan/3 = 30, Ao = and 
\x > 0. The corresponding contour plots of BF(b — > sf), and Q.h 2 are shown in 
right frame of the figure. One can see that most of the area displayed is excluded. In 
this case, slepton masses are quite light, in the vicinity of a few hundred GeV. As m 1/ / 2 
increases, ultimately m^ becomes greater than m^, and one violates the cosmological 
constraint on stable charged relics from the Big Bang. On the other hand, of the surviving 
parameter space, there are large regions with relatively low % 2 . The plot with mo(l) = 
100 GeV has a rather broad band of low % 2 . In this case, neutralinos in the early 
universe can annihilate by a combination of ^-channel slepton exchange (as in the bulk 



region of mSUGRA), and by neutralino-slepton co-annihilation. In addition, smuons 
and mu sneutrinos are relatively light, giving a large, positive contribution to Aa^, 
while top squarks inhabit the TeV and beyond range, effectively suppressing anomalous 
contributions to BF(b — > sy) . In addition to the region with low values of mo(3) and low 
mi/2, an important portion the HB/FP region survives since low mo(l) value provides 
big enough contribution to Aa^ while big mo(3) value keeps SUSY contribution to 
BR(b — > sy) suppressed. 

The scenario of NMH SUGRA 
model works if mo(l) — mo (2) <C 
mo (3), and leads to spectra typically 
with squarks and third generation 
sleptons in the TeV range, while first 
and second generation sleptons have 
masses in the range of 100-300 GeV. 
The presence of rather light first and 
second generation sleptons in the 
sparticle mass spectrum in general 
leads to enhancements in leptonic 
cross sections from superparticle 
production at collider experiments, 
compared to the case where selec- 
trons and smuons are in the multi- 
TeV range — sleptons may now be 
produced with non-negligible cross 
sections at colliders, and also in that 
their presence enhances the leptonic 
branching fractions of charginos and 
especially neutralinos. Therefore it 
is worth to discuss the implications 
of light selectrons and smuons for 
the Fermilab Tevatron collider, the 
CERN LHC and a linear e + e~ collides^ 

It is usually expected that pp — ► W*W{~X and W\Z^X will be the dominant production 

cross sections at the Tevatron if sparticle are accessible. If W\ — > £V(Z\ and Z 2 — > 
HZ\ , then clean trilepton signals may occur at an obser vable rate. Signal and background 
rates have recently been investigated in Ref. lll2lLll22Lll23h . where Tevatron reach plots 
may also be found. 

Fig. 1 1 presents isolated trilepton Isajet 7.69 signal after using cuts SC2 of Ref. lll22ll . 
where the backgrounds are also evaluated for for the parameter space point mo (3) = 
1400 GeV, m 1/2 = 225 GeV, A = 0, tan/3 = 10 and ji > 0. The results are plotted 
versus variation in the mo(l) parameter, the signal level needed for a 5a signal with 
10 fb 1 is also denoted. The error bars show the Monte Carlo statistical error. When 
mo(l) = mo(3) (at mo(l) = 1400 GeV), the results correspond to the mSUGRA model, 
and the isolated trilepton signal is well below discovery threshold. Atmo(l) < 200 GeV 
(which is favorable by 8a^ data) the light sleptons begin to dominate neutralino three 
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FIGURE 11. Rates for isolated trilepton events at the 
Ferm ilab Tevatron pp collider, after cuts SC2 from 

Ref. JITflfmlfEl 



body decays rates, and consequently the trilepton cross section rises steeply, to the level 
of observability. Eventually chargino and and neutralino two-body decays to sleptons 
turn on (in this case, first Z2 — > £r£), and trilepton rates become very high. 

Similarly, at the the CERN LHC, m (l) drops below about 200 GeV (the value 
typically needed to explain the (g — 2)^ anomaly), multilepton rates rise steeply. Thus, 
we would expect that SUSY - as manifested in the NMH SUGRA model - would 
be easily discovered, and what's more, the signal events would be unusually rich in 
multilepton events. Such multilepton events can be especial ly us eful for reconstructing 
sparticle masses in gluino and squark cascade decay events W124I1 . 

A linear collider operating at y/s = 0.5 — 1 TeV may be the next frontier particle 
physics accelerator beyond the CERN LHC. Depending on sparticle masses and the 
collider energy, charginos and neutralinos may or may not be accessible. However, in 
the NMH SUGRA model, light first and second generation sleptons are needed both to 
explain the (g — 2)^ anomaly, but also to enhance neutralino annihilation in the early 
universe. This means slepton masses are typically in the 100-300 GeV range, and likely 
within reach of a linear e + e~ collider llssll . 



SUSY GUTS AND YUKAWA UNIFIED SUSY MODELS 

The existence of the weak scale supersymmetry leads to gauge coupling unification at the 
scale of Mqut — 2 16 GeV which in its turn is compelling hing for SUSY grand unified 
theories or SUSY GUTs. Such a SUSY GUT theory may be the "low energy effective 
theory" that can be obtained from some more fundamental superstring theory. Models 
based on the gauge group SU(5) are compelling in that t hey expla in the apparently 
ad-hoc hypercharge assignments of the SM fermions lll25lll26lll27ll . However, many 
SU (5) SUSY GUT mode ls as formulated in four dimensions are already excluded by 
proton decay constraints [ 1 28] . SU(5) SUSY GUT models can also be formulated in 
five or more dimensions, where compactification o f the extra dimensions leads to a break 
down in gauge symmetry lll29lll30lll3lLll32Lll33h . These models can dispense with the 
unwieldy large Higgs representations required by four-dimensional models, and can also 
be constructed to suppress or eliminate proton decay entirely. 

Further step to the complete gauge and matter field family unification (as well as 
Yukawa coupling unification as we discuss below) can be achieved i n SUSY GUTs 



models b ased on SO(10) group which therefore looks very intriguing 1 134. 1135L 113 
In addition to unifying gauge couplings, 

• They unify all matter of a single generation into the 16 dimensional spinorial 
multiplet of 5O(10). 

• The 16 of SO( 10) contains in addition to all SM matter fields of a single generation 
a gauge singlet right handed neutrino state whic h naturally lead s to a mass for 
neutrinos. The well-known see-saw mechanism iflM flH fl4lll implies that if 
m Vt ~ 0.03 eV, as suggested by atmospheric neutrino data II142L [143J, then the 
mass scale associated with is very close to the GUT scale: i.e. M N ~ 10 15 GeV. 



50(10) explains the apparently fortuitous cancellation of triangle anomalies within 
the SM. 

The structure of the neutrino sector of 50(10) models lends itself to a succ essfu l 



theory of baryogenesis via intermediate scale leptogenesis (For a review, see [ 144]. 

• In the simplest 50(10) SUSY GUT models, the two Higgs doublets of the MSSM 
occupy the same 10 dimensional Higgs multiplet 0(10). The superpotential then 
contains the term / 9 f $(16)$ (16) (j) (10) H where / is the single Yukawa cou- 
pling for the third generation. Thus, the simplest 50(10) SUSY GUT models pre- 
dict t — b — % Yukawa coupling unification in addition to gauge coupling unification. 

It is possible to calculate the t, b and T Yukawa couplings at Q = m wea ^, and extrap- 
olate them to Mquj in much the same way one checks the theory for gauge coupling 
unification. Yukawa coupling unification turns out to depend on the entire spectrum 
of sparticle masses and mixings since these enter into the weak scale supersymmetric 
threshold corrections. Thus, the requirement of t — b — % Yukawa coupling unification 
can be a pow erful constraint on the soft SUSY breaking terms of the low energy effective 
theory IM^ fT^ fT^ fl4^ fT^ fT5fill . 

In Ref. ( 115 lfl . it was found that models with Yukawa coupling unification with R < 
1.05 (good to 5%) could be found if additional D-term splittings of scalar masses were 
included. Here, R = max(f t , f x )/min(f t , f b , f T ), where all Yukawa couplings are 
evaluated at the GUT scale. The D-term splittings occur n atura lly when the 50(10) 
gauge symmetry breaks to SU (5), and they are given by (see [1 52] and references in it). 

m 2 Q = m 2 E =ml= m\ 6 +M^ 

m D = m L = m i6 - 3M £» 
m 2 N = m 2 6 + 5Mj h 

m 2 Hud =m 2 l0 T2M 2 , (18) 

where Mq parameterizes the magnitude of the D-terms. Owing to our ignorance of the 
gauge symmetry breaking mechanism, Afp can be taken as a free parameter, with either 
positive or negative values. \Md\ is expected to be of order the weak scale. Thus, the 
D-term (DT) model is characterized by the following free parameters, 

m 16 , mio, Ml, m 1/2 , A , tan/3, sign(ju). (19) 

Using the DT model, Yukawa unification good to 5% was found when soft term 



parameters m\(, and m\Q were scanned up to 1.5 TeV, but only for ji < values [151, 
The essential quality of the D-term mass splitting is that it gave the value of m# t( a head 
start over m# rf in running towards negative values, as is required for REWSB. Good 
relic density was also found in the stau co-annihilation and hyperbolic branch/focus 
point (HB/FP) region, but at some cost to the degree of Yukawa coupling unification. 

In Ref. 01531], it was found that Yukawa coupling unification good to only 30% could 
be achieved in DT models with /i > when mi6 values were scanned up to 2 TeV. The 
models with the best Yukawa coupling unification were found to have soft term relations 



A 2 Q ~ 2m 2 w ~ 4mf 6 , 



(20) 



which had also been found by Bagger et al. in t he context of r adiatively driven inverted 
scalar mas s hierarchy (IMH) models lfl5l EH EH EH EH EH • 

InRef. [EH, a model with just GUT scale Higgs mass splittings was also examined 
(HS), while all other scalars remained universal. 

The parameter space of the HS model is that of Eq. [T9J but where the D-term splitting 
is only applied to the last of the relations in Eq. [TSl Yukawa coupling unification in the 
HS model was found to be c omp arable to the DT model case when values up to 
2 TeV were scanned. In Ref. (Auto et al), soft term values of mid U P t0 20 TeV 
were explored. In this case, Yukawa unified solutions to better than 5% were found for 
/i > for the HS model when very large values of m\(, > 5 — 10 TeV were scanned. 
The large scalar masses that gave rise to Yukawa unification also acted to suppress 
neutralino annihilation in the early universe, so that rather large values of the relic 
density were found. Models with f^/* 2 < 0.2 could be found, but only at the expense 
of accepting Yukawa coupling unification to 20%, rather than 5%. The models found 
with low relic density generally had either a low \x value, or were in the light Higgs 
annihilation corridor, with 2mg [ ~ w%. 



Similar work has been carried on by Blazek, Dermisek and Raby (BDR). In Ref. [161 



1 16211 . the BDR group used a top-down approach to the RG sparticle mass solution to find 
Yukawa unified solutions for \x > 0, where they also noted that in this case the HS model 
worked better than the DT model. In their approach, the third generation fermion masses 
and other electroweak observables were an output of the program, so that starting with 
models with perfect Yukawa coupling unification, they would look for solutions with a 
low x 2 value constructed from the low energy observables. The BDR Yukawa unified 
solutions were also characterized by soft term IMH model boundary conditions. The 
solutions differed from those of Ref. lll60ll in that they always gave a very low value of 
Higgs mass yua and also small \x parameter, indicative of a mixed higgsino-bino LSR 
In Ref. W163I1 . the neutralino relic density was examined for the BDR solutions. Their 
low jU and ma values generally led to very low values of un l ess m i/2 was small 

enough compared to \x that the LSP was in th e mixed higgsino-bino region. 

Finally, one more step was done in to solve the problem of tension between 

Yukawa unification and low relic density in SO(10) SUSY GUTs models. In this paper, 
two methods were explored to reconcile Yukawa unified sparticle mass solutions with 
the neutralino relic density. The first case is to allow splitting of the third generation 
of scalars from the first two generations. By decreasing the first and second generation 
scalar masses, sparticle mass solutions with very light ur and cr squark masses, in the 
90-120 GeV range were obtained. The other examined scenario was the case of non- 
degenerate gaugino masses at the GUT scale. Beginning with any of the Yukawa unified 
solutions found in an earlier study, it was found that by dialing M\ to large enough 
values, a (partially) wino-like LSP can be generated, with a relic density in accord with 
WMAP allowed values. 

In Fig. 12, we show the example of such solution for the first scenario in the 
mi6(l) vs. m\ji plane for fixed mig(3) = 7826.5 GeV. The blue region at low m\ji 
is excluded by the LEP2 bound on chargino mass > 103.5 GeV. The black region 

on the left is excluded because the m 2 ^ squared mass is driven tachyonic, resulting in 
a color symmetry violating ground state. The green region denotes sparticle mass spec- 



trum solution with relic density 0.094 < H^/z 2 < 0.129, within the WMAP favored 
regime, while the yellow region denotes even lower relic density values, wherein the 
CDM in the universe might be a mixture of neutralinos plus some other species. The red 
region has % h 2 < 0.5. The remaining unshaded regions give sparticle spectra solutions 

with ^2^ 2 > 0-5, and would be excluded by WMAP. We also show several contours of 

BF(b — > sf), Aa s ^ SY and R, the max/min ratio of GUT scale Yukawa couplings. We 

see that all constraints on the sparticle mass spectrum are within allowable limits for the 

shaded region to the right of the excluded region. 

The origin of the light squarks comes 

hsso(w) model: m 16 (3)=7826.5Gev, m 1d =9652.3Gev, M,/m 1s (3)=o.36988 ^ ^ the S term in the one-loop RGEs, 
ta^5i.ie5,A < /m 16 (3)=-2.i2432, v .>o,m t =i8oa e v which is non-zero and large in the case 

of multi-TeV valued split Higgs masses. 

A search for just two light squarks at the 

Fermilab Tevatron collider using the new 

Run 2 data should be able to either verify or 

disprove this scenario. In addition, the light 

squarks give rise to large rates for direct 

detection of dark matter, and should give 

observable rates for searches at CDMS II. 

CONCLUSIONS 

e for leaving several other interesting SUSY 
iew. Among them are CP and flavor violating 
SUS^ ph,ysi£fe S,y^Y,.R-pifity. 2 YioJ.aJion 24 beyond the MSSM (nMSSM) scenarios and 
various SUSY-rela^ed^^s^g^s of baryogenesis. There are also open fundamental prob- 
lems, which still need to be solved. One of the most important problems is Cosmological 
ESSiiKfi pV%bleMlawdAl®apwfetemiOi&©rMin of the /i term. 

^^rrR^f^fW^M^ffi flK$feffl& was paid to SUGRA models which are very 
generational noa-umyersality from 152L T We_ . ^,t^» . • i 

Wmmmm^^^m^pHW^ parameter space CDM constraints play a 

grj^iaj^ole ^gagii^^n^a/^e^^estrjc^ed regions to be tested by other experiments. 

^f)%{gi— ^Xtcaaajifaintsireatfe focus point, funnel and stau-co-annihilation 

)n8§u^f^elhelJiCtc^M#rr^antgiYqtiely cover the funnel region and (almost) all 





excluded by. LJ 
fffi^e WU _ 

grM^pompkmftBJflr^f rxi^/zef xfeeflJoaiKbtfidirect DM search experiments which has the 
p^mtklftfl,ahw£ffeQamptetel^foicaMfpja leaving 'no escape' for mSUGRA with 

ft^MoHiing experiments combined! 

There could be an indication coming from from 8a^ data that normal mass hierarchy 
in SUGRA scenario could be preferred by nature and that, may be, mSUGRA scenario 
is too simple to be true. 

Talking about GUTs models, SO(10) SUSY GUTs looks very attractive from both, 
theoretical and phenomenological points of view, and predicts very specific particle 
spectra testable experimentally. 

Finally, I would like to express my belief that the our exciting era of upcoming 



experiments will be indeed successful in hunting for Supersymmetry which could be 
just around corner! 
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